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Density of states, dynamic (optical) conductivity and phase diagram of strongly correlated and 
strongly disordered paramagnetic Anderson-Hubbard model are analyzed within the generalized dy- 
namical mean field theory (DMFT+E approximation). Strong correlations are accounted by DMFT, 
while disorder is taken into account via the appropriate generalization of self-consistent theory of 
localization. The DMFT effective single impurity problem is solved by numerical renormalization 
group (NRG) and we consider the three-dimensional system with semi-elliptic density of states. 
Correlated metal, Mott insulator and correlated Anderson insulator phases are identified via the 
evolution of density of states and dynamic conductivity, demonstrating both Mott-Hubbard and 
Anderson metal-insulator transition and allowing the construction of complete zero-temperature 
phase diagram of Anderson-Hubbard model. Rather unusual is the possibility of disorder induced 
Mott insulator to metal transition. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.30.+h 
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I. INTRODUCTION 

The importance of the electronic interaction and ran- 
domness for the properties of condensed matter is well 
knowni. Coulomb correlations and disorder are both 
driving forces of metal-insulator transitions (MITs) con- 
nected with the localization and delocalization of parti- 
cles. In particular, the Mott-Hubbard MIT is caused by 
the electronic repulsion 2 , while Anderson MIT is due to 
random scattering of non-interacting particles 3 -. Actu- 
ally, disorder and interaction effects are known to com- 
pete in many subtle ways 1 - 4 - and this problem becomes 
much more complicated in the case of strong electron 
correlations and strong disorder, determining the physi- 
cal mechanisms of Mott- Anderson MIT 1 . 

The cornerstone of the modern theory of strongly 
correlated systems is the dynamical mean-field theory 
(DMFT)££££, constituting a non-perturbative theoret- 
ical framework for the investigation of correlated lattice 
electrons with a local interaction. In this approach the 
effect of local disorder can be taken into account through 
the standard average density of states (DOS) 9 - in the 
absence of interactions, leading to the well known co- 
herent potential approximation^, which does not de- 
scribe the physics of Anderson localization. To over- 
come this deficiency Dobrosavljevic and Kotliar— for- 
mulated a variant of the DMFT where the geometrically 
averaged local DOS was computed from the solutions of 
the self-consistent stochastic DMFT equations. Subse- 
quently, Dobrosavljevic et al*£ incorporated the geomet- 
rically averaged local DOS into the self-consistency cycle 
and derived a mean-field theory of Anderson localization 
which reproduced many of the expected features of the 
disorder-driven MIT for non-interacting electrons. This 
approach was extended by Byczuk et a/. 13 to include 
Hubbard correlations via DMFT, which lead to highly 
non trivial phase diagram of Anderson-Hubbard model 
with correlated metal, Mott insulator and correlated An- 
derson insulator phases. The main deficiency of these ap- 



proaches, however, is inability to calculate directly mea- 
surable physical properties, such as conductivity, which 
is of major importance and defines MIT itself. 

At the same time the well developed approach of self- 
consistent theory of Anderson localization, based on solv- 
ing the equations for the generalized diffusion coefficient, 
demonstrated its efficiency in non-interacting case long 
ago 14 i 15 i 16 ' 17 i 18 ' 19 and some attempts to include interac- 
tion effects within this approach were also undertaken 
with some promising result o 17 i 20 . However, up to now 
there were no attempts to incorporate this approach to 
the modern theory of strongly correlated electronic sys- 
tems. Here we undertake such research, studying both 
Mott-Hubbard and Anderson MITs via direct calcula- 
tions of both the average DOS and dynamic (optical) 
conductivity. 

Our approach is based on recently proposed gener- 
alized DMFT+E approximation 21 - 2 ^ 2 ^, which on the 
one hand retains the single-impurity description of the 
DMFT, with a proper account for local Hubbard-like 
correlations and the possibility to use impurity solvers 
like NR G 25 ' 26 ! 27 , while on the other hand, allows one to 
include additional (either local or non-local) interactions 
(fluctuations) on a non-perturbative model basis. 

Within this approach we have already studied both 
single- and two-particle properties, of two-dimensional 
Hubbard model, concentrating mainly on the problem 
of pseudogap formation in the density of states of the 
quasiparticle band in both correlated metals and doped 
Mott insulators, with an application to superconducting 
cuprates. We analyzed the evolution of non-Fermi liquid 
like spectral density and ARPES spectra 2 ^, "destruction" 
of Fermi surfaces and formation of Fermi "arcs"—, as well 
as pseudogap anomalies of optical conductivity 24 . Briefly 
we also considered impurity scattering effects 23 . 

In this paper we apply our DMFT+E approach for cal- 
culations of the density of states, dynamic conductivity 
and phase diagram of strongly correlated and strongly 
disordered three-dimensional paramagnetic Anderson- 
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Hubbard model. Strong correlations are again accounted 
by DMFT, while disorder is taken into account via the 
appropriate generalization of self-consistent theory of lo- 
calization. 

The paper is organized as follows: In section [IT] we 
present a short description of our generalized DMFT+E 
approximation with application to disordered Hubbard 
model. In section Hill we present basic DMFT+E expres- 
sions for dynamic (optical) conductivity and formulate 
appropriate self-consistent equations for the generalized 
diffusion coefficient. Computational details and results 
for density of states and dynamic conductivity are given 
in section HVl where we also analyze the phase diagram of 
strongly disordered Hubbard model, following from our 
approach. The paper is ended with a short summary sec- 
tion [V] including a discussion of some related problems. 



II. BASICS OF DMFT+S APPROACH 

Our aim is to consider non-magnetic disordered 
Anderson-Hubbard model (mainly) at half-filling for 
arbitrary interaction and disorder strengths. Mott- 
Hubbard and Anderson MITs will be investigated on an 
equal footing. The Hamiltonian of the model under study 
is written as: 



H 



^eiriia + U^2n^nn, (1) 



where t > is the amplitude for hopping between near- 
est neighbors, U is the on-site repulsion, ni a — al a a ia is 
the local electron number operator, (aj CT ) is the an- 
nihilation (creation) operator of an electron with spin a, 
and the local ionic energies ej at different lattice sites are 
considered to be independent random variables. To sim- 
plify diagrammatics in the following we assume Gaussian 
probability distribution for ef. 
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Here the parameter A is just a measure of disorder 
strength, and Gaussian ("white" noise) random field of 
energy level at lattice cites produces "impurity" scat- 
tering, leading to the standard diagram technique for cal- 
culation on the averaged Green's functions^. 

DMDF+E approach was initially propose d 21 ! 22 ' 23 as 
a simple method to include non-local fluctuations, es- 
sentially of arbitrary nature, to the standard DMFT. In 
fact it can be used to include into DMFT any additional 
interaction in the following way. Working at finite tem- 
peratures T we write down Matsubara "time" Fourier 
transformed single-particle Green function of the Hub- 
bard model as: 



G(ie,p) 
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ie + fi — e(p) - S(ie) - S p («e„ 



where e(p) is the single particle spectrum, corresponding 
to free part of fl]), /i is chemical potential fixed by elec- 
tron concentration, and E(ie) is the local contribution to 
self-energy due to Hubbard interaction, of DMFT type 
(surviving in the limit of spatial dimensionality d — ► oo), 
while E p (ie) is some additional (in general momentum 
dependent) self-energy part. This last contribution can 
be due e.g. to electron interactions with some "addi- 
tional" collective modes or order parameter fluctuations 
within the Hubbard model itself. But actually it can 
be due to any other interactions (fluctuations) outside 
the standard Hubbard model, e.g. due to phonons or 
to random impurity scattering, when it is in fact local 
(momentum independent). The last interaction will be 
our main interest in the present paper. Basic assumption 
here is the neglect of all interference processes of the local 
Hubbard interaction and "external" contributions due to 
these additional scatterings (non-crossing approximation 
for appropriate diagrams) 22 , as illustrated by diagrams 
in Fig. EJ 

The self-consistency equations of generalized 
DMFT+E approach are formulated as follow s 21 ' 22 : 

1. Start with some initial guess of local self-energy 
E(ie), e.g. E(ie) = 0. 

2. Construct E p (ze) within some (approximate) 
scheme, taking into account interactions with "ex- 
ternal" interaction (impurity scattering in our case) 
which in general can depend on E(iw) and fi. 

3. Calculate the local Green function 
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N ^ie + fi- e(p) - S(ie) - E p (ie) ' 

4. Define the "Weiss field" 

S -» = E(ze) + G^». 



(4) 



(5) 



5. Using some "impurity solver" calculate the single- 
particle Green function Gd(is) for the effective An- 
derson impurity problem, placed at lattice site i 
and defined by effective action which is written, in 
obvious notations, as: 

S e g = - / dri / dr 2 c lo -(ri)^ 1 (Ti - r 2 )c+ (r 2 ) 
Jo Jo 

+ / dTUn^(T)n zl (T). (6) 
Jo 

Actually, in the following, for the "impurity solver" 
we use NR G 25 ' 26 ' 27 , which allows us to deal also 
with real frequencies, thus avoiding the complicated 
problem of analytical continuation from Matsubara 
frequencies. 



-,£ = 7rT(2n+l), 
(3) 



6. Define a new local self-energy 

E(iw) = {ioj) — G^^uj). 



(7) 
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7. Using this self-energy as "initial" one in step 1, 
continue the procedure until (and if) convergency 
is reached to obtain 

G u (ie)=G d (ie)- (8) 

Eventually, we get the desired Green function in the form 
of |3]), where E(ze) and E p (je) are those appearing at the 
end of our iteration procedure. 

For E p (ie„) of the random impurity problem we shall 
use the simplest possible one-loop contribution, given by 
the third diagram of Fig. [1] (a), neglecting "crossing" 
diagrams like the fourth one in Fig. [1] (a), i.e. just the 
self-consistent Born approximation^ 9 -, which in the case 
of Gaussian disorder @ leads to the usual expression: 

E p (ie) = A 2 G ( l£ , P) = Zim P (ie) (9) 
p 

which is actually p-independent (local). 



III. DYNAMIC CONDUCTIVITY IN DMFT+E 
APPROACH 

A. Basic expressions for optical conductivity 

Physically it is clear that calculations of dynamic con- 
ductivity are the most direct way to study MITs, as its 
frequency dependence along with static value at zero fre- 
quency of an external field allows the clear distinction 
between metallic and insulating phases (at zero temper- 
ature T = 0). 

To calculate dynamic conductivity we use the general 
expression relating it to retarded density-density corre- 
lation function q) 14 ' 19 : 

aH = -hm^ X >,q), (10) 

where e is electronic charge. 

Next we briefly follow the derivation presented in de- 
tail in Ref.— for the pseudogap problem, with necessary 
modifications for the present case. Consider full polar- 
ization loop graph in Matsubara representation shown 
in Fig. [21 which is conveniently (with explicit frequency 
summation) written as: 

q) = ( iuJ > «l) = £ *te(«*>> q) (11) 

ee' e 

and contains all possible interactions of our model, de- 
scribed by the full shaded vertex part. Actually we im- 
plicitly assume here that simple loop contribution with- 
out vertex corrections is also included in Fig. [2j which 
shortens further diagrammatic expressions 2 — Retarded 
density-density correlation function is determined by ap- 
propriate analytic continuation of this loop and can be 



written as**: 

X%;,q) = J ^ — {[/( e+ ) - f(e-)] $f A (q,^ 

+f(e-)*™(q 1 w) - /( £+ )0>^(q, W )} , (12) 

where /(e) - Fermi distribution, e± = e ± while 
two-particle loops $f A (q,w), $f fl (q, oj), ^^(q,^) are 
determined by appropriate analytic continuations (ie + 
iu — > e + oj + i5,ie — > e ± i5,S — > +0) in (fTTj) . Then we 
can write dynamic (optical) conductivity as: 

- /( £ _)] [$? A (q,oj) $? A (0,u)} 
+/(e_)[$^(q, w )-$f>,c)] 
-/(e + )[$^(q,c)-<i>^(0,c)]}, (13) 

where the total contribution of additional terms with zero 
q can be shown (with the use of general Ward identities) 
to be zero. 

To calculate $i e j e ' (iuj, q), entering the sum over Mat- 
subara frequencies in pip , in DMFT+E approxima- 
tion, which neglects interference between local Hub- 
bard interaction and impurity scatterings, we can write 
down Bethe-Salpeter equation, shown diagrammatically 
in Fig. [3l where we introduce irreducible (local) vertex 
Ui £ i £ > (iuj) of DMFT and "rectangular" vertex, containing 
all interactions with impurities. Analytically this equa- 
tion can be written as: 

+<P Q iE (iuj, q) Uieie" {iuj)<5>. l£ n ie , (iu, q), (14) 

e" 

where $° e (iw,q) is the desired function calculated ne- 
glecting vertex corrections due to Hubbard interaction 
(but taking into account all interactions due to impurity 
scattering). Note that all g-dependence here is deter- 
mined by &® e (iuj, q) as the vertex Ui E i S '(iuj) is local and 
q- independent. 

As we noted in Refill, it is clear from (fT3")) that cal- 
culation of conductivity requires only the knowledge of 
g 2 -contribution to $(ia;,q). This can be easily found in 
the following way. First of all, note that all the loops 
in (|14p contain g-dependence starting from terms of the 
order of q 2 . Then we can take an arbitrary loop (crossec- 
tion) in the expansion of (|14[) (see Fig. [3]), calculating 
it up to terms of the order of q 2 , and make resumma- 
tion of all contributions to the right and to the left from 
this crossection, putting q = in all these graphs. This 
is equivalent to simple ^-differentiation of the expanded 
version of Eq. (fT4")) . This procedure immediately leads 
to the following relation for ^-contribution to (fTTj) : 

= km — '—^ '- = 2^ T&(tw, q = (M e (iw) 

<j— >u q- * — ' 

(15) 
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where 

&(fao = nm m^m^s (16) 

where $? e (iu;,q) contains vertex corrections only due to 
impurity scattering, while one-particle Green's functions 
entering it are taken with self-energies due to both im- 
purity scattering and local DMFT-like interaction, like 
in Eq. ((3J. The vertex ^i e {iw, q = 0) is determined dia- 
grammatically as shown in Fig. 31 or analytically: 

J iE (iuj, q = 0) = 1 + ^ Uieie"(iu)® iE ><ie>(iu,<l = 0). 
e' e" 

(17) 

Now using Bethe-Salpeter equation (fT4"]) we can write 
explicitly: 

7 fe (wj,q = 0) = 1 + 

$ i£fe / (iw, q = 0) - $^(i£J, q = 0) 



^ $° e (ia,,q=0) 
$? E (iw,q = 0) 



(18) 



For q = we have the following Ward identity, which can 
be obtained by direct generalization of the proof given 



,14,28 



(see details in Appendix of Rcf. - 



(-iu>)$ie(iuj, q = 0) = (-iu) ^ ^ iei£ ' (iw, q = 0) 

e' 

= 5^G(ie + iw,p)-5^G(ie J p). (19) 



Denominator of (|18|1 contains vertex corrections only 
from impurity scattering, while Green's functions here 
are "dressed" both by impurities and local (DMFT) Hub- 
bard interaction. Thus we may consider the loop entering 
the denominator as dressed by impurities only, but with 
"bare" Green's functions: 



G (ie,p) 



1 



ie + n — e(p) — E(ze) ' 



(20) 



where E(ie) is local contribution to self-energy from 
DMFT. For this problem we have the following Ward 
identity, similar to ([13)) (see Appendix of Ref. 24 ): 

G ( ie + iuJ i P) - G ( i£ > P) = 
p p 

= $i e (iuj, q = 0) [E(ie + iu) - £(ie) - iu] = 

= $?(»w,q=0)[AS(« I ;)-»w], (21) 



where we have introduced 

AE(iw) = S(i£ + iw) - E(is). 



(22) 



Thus, using (JTjJJ), (|2T|) in (|T8]) we get the final expression 
for 7 i£ (icj,q = 0): 



7ie(iw,q = 0) = 1 



AE(iw) 



(23) 



Then (|15)) reduces to: 

0M = VCH 



i 



AE(iw) 



?,u; 



(24) 



Analytic continuation to real frequencies is obvious and 
using (fT5)) , in |T3]) we can write the final expression 
for the real part of dynamic (optical) conductivity as: 

ReaH = ^ f de[f{e-) - f(e + )} 




E*( E+ )-E^( £ _) 



£*( £+ )-E«( £ _) 



(25) 



Thus we have achieved a great simplification of our prob- 
lem. To calculate dynamic conductivity in DMFT+E ap- 
proximation we only have to solve single-particle prob- 
lem as described by DMFT+E procedure above to deter- 
mine self-consistent values of local self-energies E(e±), 
while non-trivial contribution of impurity scattering are 
to be included via (fT6|) . which is to be calculated in 
some approximation, taking into account only interac- 
tion with impurities (random scattering), but using the 
"bare" Green's functions of the form (|2"0)) , which in- 
clude local self-energies already determined via the gen- 
eral DMFT+E procedure. Actually ([25]) provides also 
an effective algorithm to calculate dynamic conductivity 
in standard DMFT (neglecting impurity scattering), as 
(I16p is then easily calculated from a simple loop diagram, 
determined by two Green's functions and free scalar ver- 
tices. As usual, there is no need to calculate vertex cor- 
rections within DMFT itself, as was proven first consid- 
ering the loop with vector vertices^. Obviously, Eq. 
(|25|) provides effective interpolation between the case of 
strong correlations without disorder and the case of pure 
disorder, without Hubbard correlations, which is of major 
interest to us. In the following we shall see that calcula- 
tions based on Eq. ([25]) give a reasonable overall picture 
of MIT in Anderson-Hubbard model. 



B. Self-consistent equations for generalized 
diffusion coefficient and conductivity 

Now to calculate optical conductivity we need the 
knowledge of the basic block $° e («tc, q), entering (fl"6|) . 
or, more precisely, appropriate functions analytically con- 
tinued to real frequencies: $® ra (lj, q) and &° RR (uj,q), 
which in turn define (f>® RA (uj) and (jP e RR {ui) entering (|25p . 
and defined by obvious relations similar to (|16j) : 



iORA 



M = lim , 



4> 



ORA 



(cj,q)-^ RA (u,0) 



lim 

q^0 



$f fl (w,q)-$™( W ,0) 



(26) 



(27) 
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By definition we have: 



xr^( £ _,p_;e + ,p + ) (28) 



*°* fl ( W) q) =^G fl ( £+ , P+ )G fi (e_,p_) 



xr M (e_,p_; £+ ,p + ), (29) 



which are shown diagrammatically in Fig. [5] and 
p± = p ± #. Here Green's functions G fl (£ + ,p+) and 
G j4 (e_,p_) are defined by analytic continuation (is — > 
e ± i<5) of Matsubara Green's functions ([3]) determined 
via our DMFT+E algorithm (gj) - ©, while vertices 
T RA (e^, p ; e + , p + ) and r flfl (e_, p ; e + , p + ) contain all 
vertex corrections due to impurity scatterings. 

The most important block <&® RA (u>, q) can be calcu- 
lated using the basic approach of self-consistent theory 
of localizatio n 14 ' 15 i 16 ' 17 i 18 ' 19 with appropriate extensions, 
taking into account the role of the local Hubbard inter- 
action using DMFT+S approach. The only important 
difference with the standard approach is that equations 
of self-consistent theory are now derived using 



G R > A (e,p) 



1 



£ + )U - e (p)-S^( e )-Sf4( £ ) 



(30) 



containing DMFT contributions T, R ' A (e), not only impu- 
rity scattering contained in: 

Ef m >) = A 2 G R,A {e, p) = ReSWe) ± i 7 ( e ) (31) 



where 7(e) = 7rA 2 A^(e) and N(e) is the density of 
states renormalized by Hubbard interaction, accounted 
via DMFT+E and given as usual by: 



TT * J 



(32) 



Following all the usual steps of standard 
derivatio n 14 ' 15 i 16 ' 17 i 18 ' 19 we obtain diffusion like (at 
small ui and q) contribution to $>® RA (u>, q) as: 



^(q,a>) 



2iriN(e) 
Co + iD(uj)q 2 



(33) 



where important difference with the single-particle case 
is contained in 

£ = e+-e_ -£ fl (e+) + £ A (e_) = 
= w-E Ji (£ + ) + S i ( £ _)=^-AS RA H (34) 

which replaces the usual u) term in the denominator of the 
standard expression for <&® RA (u>, q). On general grounds 



it is clear that in metallic phase for u) — > we have 
AS M (w = 0) = 2ilm£(e) ~ Max{T 2 ,e 2 }, reflecting 
Fermi-liquid behavior of DMFT (conserved by elastic im- 
purity scattering). At finite T it leads to the usual phase 
decoherence due to electron - electron scattering^ 4 -. Gen- 
eralized diffusion coefficient D(lu) will be determined by 
solving the basic self-consistency equation, introduced 
below. 

Now using (|33"1) in (|26[) we easily obtain: 



2kN{e)D{uj) 



uj 2 11- 



Ag« A (^) 



(35) 



Then using (J35J) in (S5J) , for ui -> and for T = we get 
just the usual Einstein relation for static conductivity: 



cr(0) = e 2 N(0)D(0) 



(36) 



All contributions form Hubbard interaction are reduced 
to renormalization of the density of states at the Fermi 
level and also of diffusion coefficient D(0). 
Then reduces to: 



ReaH = ^ J°° de [/(e_) - /(e+)] x 



Re 



AY, rr (lo) 



(37) 



where the second term actually can be neglected at small 
ui, or just calculated from (j2"T|) taking &® RR (cu,q) given 
by the usual "ladder" approximation (|A10|) . 

Now we have to formulate our basic self-consistent 
equation, determining the generalized diffusion coeffi- 
cient D(w). Again we follow all the usual steps of self- 
consistent theory of localization (see details in the Ap- 
pendix A), taking into account the form of our single- 
particle Green's function (piO)) , and not restricting analy- 
sis to small to limit. Then we can write the generalized 
diffusion coefficient as: 



D(u) 



< v 



M(uj) 



(38) 



where d is spatial dimensionality and average velocity 
< v > is defined in (|A6|) (to a good approximation it 
is just the Fermi velocity), while the relaxation kernel 
M(oj) satisfies self-consistency equation, similar to that 
derived in Refsi 14 ' 15 i 16 ' 17 i 18 ' 19 , using "maximally crossed" 
diagrams for irreducible impurity scattering vertex (built 
with Green's functions (|3"0"|0 : 



MM = -AE&LM + A 4 V(AG P ) 



1 



p 



q 



with 



AG p =G fl (e + ,p)-G A (e_,p) 



iD{w)q 2 
(39) 

(40) 
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and AEg* (w) = S^ p (e + ) - E^ ip (e_) is due to im- 
purity scattering. It is important to stress once again 
that there are no contributions to this equation due to 
vertex corrections, determined by local Hubbard interac- 
tion. Using the definition ([38]) Eq. (|39|) can be rewritten 
as self-consistent equation for the generalized diffusion 
coefficient itself: 

D(u) = i^^U - + 

p q J 

which is to be solved in conjunction with our DMFT+E 
loop Q-©. Due to the limits of diffusion approxima- 
tion, summation over q in ()41|) is to be restricted to: 

q<k = Min{r 1 ,p F } (42) 

where I =< v > /2 r y(0) is an elastic mean-free path, pp 
is the Fermi momentu m 17 ! 19 . 

Solving (|4"Tj) for different sets of parameters of our 
model and using it in (|37[) with regular contributions 
from (|A10[) we can calculate dynamic (optical) conduc- 
tivity in different phases of Anderson-Hubbard model. 

IV. RESULTS AND DISCUSSION 

We performed extensive numerical calculations for sim- 
plified version of three-dimensional Anderson-Hubbard 
model on cubic lattice with semi-elliptic DOS of "bare" 
band of the width W = 2D: 

N °( £ ) = ~^2^ D2 - £2 ( 43 ) 
■kD z 

DOS is always given in units of number of states per 
energy interval, per lattice cell volume a 3 (a is lattice 
spacing), per spin. Some related technical details are 
given in Appendix B. 

Mostly we shall concentrate on the half - filled case, 
though some results for finite dopings will be also pre- 
sented. Fermi level is always placed at zero energy. 

As "impurity solver" of DMFT we employed the re- 
liable numerically exact method of numerical renormal- 
ization group (NRG) 25 ' 26 ' 2 , Calculations were done for 
temperatures T ~ 0.001D, which effectively makes tem- 
perature effects in DOS and conductivity negligible. Dis- 
cretization parameter of NRG was always A=2, number 
of low energy states after truncation 1000, cut-off near 
Fermi energy !) -6 , broadening parameter &=0.6. 

Below we present only a fraction of most typical re- 
sults. 



A. Density of states evolution 

Within the standard DMFT approach density of states 
of the half- filled Hubbard model has a typical three peak 



structure: a narrow quasiparticle band (central peak) de- 
velops at the Fermi level, with wider upper and lower 
Hubbard bands forming at e ~ ±U/2. Quasiparticle 
band narrows further with the growth of U in metallic 
phase, vanishing at critical U C 2 ~ 1.5W, signifying the 
Mott-Hubbard MIT with a gap opening at the Fermi 
leve£&21. 

In Fig. [6] we present our DMFT+E results for the 
density of states, obtained for U — 2.5D = 1.25W typ- 
ical for correlated metal without disorder, for different 
degrees of disorder A, including strong enough values, 
actually transforming correlated metal to correlated An- 
derson insulator (see next subsection IIV Bp . As may be 
expected, we observe typical widening and damping of 
DOS by disorder. 

More unexpected are the results obtained for the val- 
ues of U typical for Mott insulator without disorder, as 
shown in Fig. for U = 4.5D = 2.25W. We see the 
restoration of central peak (quasiparticle band) in DOS 
with the growth of disorder, transforming Mott insulator 
either to correlated metal or correlated Anderson insu- 
lator. Similar behavior of DOS was recently obtained in 
Ref.— . However, in our calculations the presence of dis- 
tinct Hubbard bands was observed even for rather large 
values of disorder, with no signs of vanishing Hubbard 
structure of DOS, which was observed in of Ref^ 3 . This 
is probably due to very simple nature of our approxima- 
tion for DOS under disordering, though we must stress 
that this difference may be also due to another model of 
disorder used in Refi^ 3 -, i.e. flat distribution of e< in ([1]) 
instead of our Gaussian case ([2]). Though unimportant, 
in general, to physics of Anderson transition, the type of 
disorder may be significant for the DOS behavior. 

It is well known, that hysteresis behavior of DOS is ob- 
tained for Mott-Hubbard transition if we perform DMFT 
calculations with U decreasing from insulating phase&^i. 
Mott insulator phase survives for the values of U well 
inside the correlated metal phase, obtained with the in- 
crease of U. Metallic phase is restored at U c i w 1.014 7 . 
The values of U from the interval U c \ < U < U c i are 
usually considered as belonging to coexistence region 
of metallic and (Mott) insulating phases, with metallic 
phase being thermodynamically more stable ^ 27 ' 29 . 

In Fig. [S]we present our typical data for DOS with dif- 
ferent disorder for the same value of U = 2.5D = 1.25W 
as in Fig. but for hysteresis region, obtained by de- 
creasing U from Mott insulator phase. Again we observe 
the restoration of central peak (quasiparticle band) in 
DOS under disordering. Note also the peculiar form of 
DOS around the Fermi level during this transition - a 
narrow energy gap is conserved until it is closed by dis- 
order, and central peak is formed from two symmetric 
maxima in DOS joining into quasiparticle band. This 
reminds similar behavior observed in periodic Anderson 
model^. Apparently this effect was unnoticed in previ- 
ous calculations of DOS in coexistence region^ 7 , (in the 
absence of disorder), while in our case it was obtained 
mainly due to our use of very fine mesh of the values of 
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disorder parameter A. 

The physical reason for rather unexpected restoration 
of the central (quasiparticle) peak in DOS is in fact clear. 
The controlling parameter for its appearance or disap- 
pearance in DMFT is actually the ratio of Hubbard in- 
teraction U and the bare bandwidth W = 2D. Under 
disordering we obtain the new effective bandwidth W e f / 
(in the absence of Hubbard interaction) which grows with 
disorder, while semi-elliptic form of the DOS, with well- 
defined band edges is conserved in self-consistent Born 
approximation ©. This leads to the decrease of the ra- 
tio U/W e ff, which induces the reappearance of quasi- 
particle band in our model. This will be illustrated in 
more detail in subsection IIV C| where our DOS calcu- 
lations within DMFT+S approach for a wide range of 
parameters will be used to study the phase diagram of 
Anderson-Hubbard model. 



B. Dynamic conductivity: Mott— Hubbard and 
Anderson transitions 

Real part of dynamic (optical) conductivity was cal- 
culated for different combinations of parameters of our 
model directly from Eqs. (|37|). (|A10|) and gT]) using the 
results of DMFT+S loop © - © as an input. The val- 
ues of conductivity below are given in natural units of 
e 2 /Ha (a - lattice spacing). 

In the absence of disorder we obviously reproduce the 
results of the standard DMFT— £ with dynamic conduc- 
tivity characterized in general by the usual (metallic) 
Drude-like peak at zero frequency and wide absorption 
maximum at u> ~ U, corresponding to transitions to the 
upper Hubbard band. With the growth of U Drude peak 
decreases and vanishes at Mott transition, when only 
transitions through the Mott-Hubbard gap contribute. 
Introduction of disorder leads to qualitative changes in 
the frequency dependence of conductivity. Below we 
mainly show the results obtained for the same values of 
U and A as were used above to illustrate DOS behavior. 

In Fig. [9] we present real part of dynamic (optical) 
conductivity for half-filled Anderson-Hubbard model for 
different degrees of disorder A, and U = 2.5D, typical for 
correlated metal. Transitions to upper Hubbard band 
at uj rsj U are practically unobservable in these data. 
However, it is clearly seen that metallic Drude peak at 
zero frequency is widened and suppressed, being gradu- 
ally transformed to a peak at finite frequncies due to ef- 
fects of Anderson localization. Anderson transition takes 
place at A c w 0.7AD = 0.37VF (which in all our graphs 
(including those for DOS) corresponds to curve 3). Note 
that this value is actually dependent on the value of cut- 
off (|42p , which is defined up to a constant of the order of 
unit y 17 ! 19 . Naive expectations may have lead to a con- 
clusion, that a narrow quasiparticle band at the Fermi 
level, which forms in general case of highly correlated 
metal, may be localized much more easily than typical 
conduction band. However, we see that these expecta- 



tions are just wrong and this band is localized only at 
strong enough disorder A c ~ D, the same as for the 
whole conduction band of the width ~ W. This is in 
accordance with previous analysis of localization in two- 
band model^. 

More important is the fact that in DMFT+E approx- 
imation the value of A c is independent of U as all inter- 
action effects enter Eq. (|4ip only via AE M (w) — > for 
uj — ► (at T = 0), so that interaction just drops out at 
w = 0. This is actually the main deficiency of our ap- 
proximation, which is due to our neglect of interference 
effects between interaction and disorder scattering. Im- 
portant role of these interference effects is known for a 
long timers. However, despite the neglect of these effects 
we are able to produce physically sound interpolation be- 
tween two main limits of interest - pure Anderson transi- 
tion due to disorder and Mott-Hubbard transition due to 
strong correlations. Thus we consider it as a reasonable 
first step to future complete theory of MIT in strongly 
correlated disordered systems. 

In Fig. [TU] we present real part of dynamic (opti- 
cal) conductivity for different degrees of disorder A, and 
U = 4.5.D, typical for Mott-Hubbard insulator. At the 
insert we show our data for small frequencies, which allow 
clear distinction of different types of conductivity behav- 
ior, especially close to Anderson transition, or in Mott 
insulator phase. In this figure we clearly see the contri- 
bution of transitions to upper Hubbard band at uj ~ U. 
More importantly we observe, that the growth of dis- 
order produces finite conductivity within the frequency 
range of Mott-Hubbard gap, which correlates with the 
appearance of quasiparticle band (central peak) in DOS 
within this gap, as shown in Fig. [JJ In general case, 
this conductivity is metallic (finite in the static limit of 
uj = 0) for A < A c , while for A > A c at small frequen- 
cies we obtain Kea(uj) ~ uj 2 , which is typical of Anderson 
insulatorii^^iii^ii 9 .. Note that due to a finite inter- 
nal accuracy of NRG numerics, small but finite spurious 
contributions to IniYi R ' A (e = 0) always appear^ 7 - and for- 
mally grow with U. These contributions are practically 
irrelevant in calculations of conductivity in metallic state. 
However, in Anderson insulator these spurious terms con- 
tribute via to in Eq. (|4ip and lead to unphysical finite 
dephasing effects at uj — (or T — 0), which can sim- 
ulate small finite static conductivity. To exclude these 
spurious effects we had to make appropriate subtractions 
in our data for lmY: R ' A (e) at e = 0. 

Rather unusual is the appearance of low frequency 
peak in Reer(u;) at low frequencies even in metallic phase. 
This is due to importance of weak localization effects, as 
can be clearly seen from Fig. 1111 where we compare 
the real part of dynamic conductivity for different de- 
grees of disorder A and U — 1.5D, obtained via our 
self-consistent approach (taking into account localization 
effects via "maximally crossed" diagrams) with that ob- 
tained using the "ladder" approximation for <&® RA (oj,q) 
(similar to (|A10j) ). which neglects all localization effects. 
It is clearly seen that in this simple approximation we just 
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obtain the usual Drude-like peak at u> = 0, while the ac- 
count of localization effects produce the peak in Recr(w) 
at low (finite) frequencies. Metallic state is defined 2 
by the finite value of zero temperature conductivity at 
w = 0. 

Up to now we presented only conductivity data ob- 
tained with increase of U from metallic to (Mott) in- 
sulating phase. As we decrease U from Mott insulator- 
hysteresis of conductivity is observed in coexistence re- 
gion, defined (in the absence of disorder, A = 0) by 
U c i < U < U C 2- Typical data are shown in Fig. [T2J 
where we present the real part of dynamic conductiv- 
ity for different degrees of disorder A and U = 2.5D, 
obtained from Mott insulator with decreasing U, which 
should be compared with those shown in Fig. [9] Tran- 
sition to metallic state via the closure of a narrow gap, 
"inside" much wider Mott-Hubbard gap, is clearly seen, 
which correlates with DOS data shown in Fig. [5] 



C. Phase diagram of half— filled Anderson—Hubbard 
model. 

Phase diagram of half-filled Anderson-Hubbard model 
was studied in Ref^, using the approach, based on di- 
rect DMFT calculations for a set of random realizations 
of site energies £j in ([1]) with subsequent averaging to 
get both the standard average DOS and also geometri- 
cally averaged local DOS, which was used to determine 
transition to Anderson insulator phase. Here we present 
our results for the zero-temperature phase diagram of 
half-filled paramagnetic Anderson-Hubbard model, ob- 
tained from extensive calculations of both average DOS 
and dynamic (optical) conductivity in DMFT+S approx- 
imation. It should be noted, that conductivity calcula- 
tions are the most direct way to distinguish metallic and 
insulating phases^. 

Our phase diagram in the disorder-correlation (A, U)— 
plane is shown in Fig. 1131 Anderson transition line A c ss 
0.37W = 0.7AD was determined as the value of disorder 
for which the static conductivity becomes zero at T = 0. 
Mott-Hubbard transition may be determined either via 
the disappearance of central peak (quasiparticle band) in 
DOS, or from conductivity, e.g. looking for the closure of 
gap in dynamic conductivity in insulating phase, or from 
vanishing static conductivity in metallic region. All these 
methods were used and appropriate results are shown for 
comparison in Fig. [13] 

We already stressed that DMFT+E approximation 
gives the universal ([/-independent) value of A c . This is 
due to neglect of interference between disorder scattering 
and Hubbard interaction, and leads to the main (over) 
simplification of our phase diagram, compared with that 
obtained in Refi^. Note, that direct comparison of our 
critical disorder value with those of Ref is complicated 
by different types of random site-energies distributions 
used here (Gaussian) and in Ref^ (rectangular distri- 
bution). As a rule of thumb (cf. second reference ir>i£) 



our Gaussian value of A c should be multiplied by yl2 to 
obtain the critical disorder value for rectangular distri- 
bution. This gives A c w 1.28 in rather good agreement 
with A C (U = 0) w 1.35W value of RefJ^, justifying our 
cutoff choice in (|42|) . 

The influence of disorder scattering on Mott-Hubbard 
transition is highly non-trivial and in some respects is in 
qualitative agreement with results of RefJ^. The main 
difference is that our data indicate the survival of Hub- 
bard bands structures in DOS even in the limit of rather 
large disorder, while these were claimed to disappear in 
Refi^. Also we obtain the coexistence region smoothly 
widening with the growth of disorder and not disappear- 
ing at some "critical" point, as in Ref. 13 . The borders 
of our coexistence region, which in fact define the bor- 
ders of Mott insulator phase obtained with increasing or 
decreasing U, are determined by the lines of U c i(A) and 
C/ C 2(A) shown in Fig. [T3] which are obtained from the 
simple equation: 

Ucl,c2(A) _ Ucl,c2 

W eff W 1 ' 

with 

W eff = W^fl + W^ (45) 

which is the effective bandwidth in the presence of disor- 
der, calculated for U = in self-consistent Born approx- 
imation ([9]). Thus the borders of coexistence region are 
given by: 



C/ cl , c2 (A) = U cl , c2 \j 1 + 16^2 (46) 

which are explicitly shown in Fig. EH by dotted and full 
lines, defining the borders of Mott insulator phase. Nu- 
merical results for disappearance of quasiparticle band 
(central peak) in DOS, as well as points following from 
qualitative change in conductivity behavior, are shown 
in Fig. [I3]by different symbols demonstrating very good 
agreement with these lines, confirming the ratio (J44J) as 
controlling parameter of Mott transition in the presence 
of disorder. 

Most striking result of our analysis (also qualitatively 
demonstrated in Refji 3 -) is the possibility of metallic state 
being restored from Mott-Hubbard insulator with the 
growth of disorder. This is clear from the phase diagram 
and is nicely demonstrated by our data for (static) con- 
ductivity shown in Fig. [14] for several values of U > U C 2 
and disorder values A < A c . At the insert in Fig. [T4"l we 
also illustrate static conductivity hysteresis, observed in 
coexistence region of the phase diagram, obtained with 
U decreasing from Mott insulator phase. 

D. Doped Mott insulator 

All results presented above were obtained for half- 
filled case. Here we briefly consider deviations from 
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half-filling. In metallic phase, doping from half-filling 
does not produce any qualitative changes in conductiv- 
ity behavior, which only demonstrates Anderson transi- 
tion with the growth of disorder. Thus we shall concen- 
trate only on the case of doped Mott insulator. Strictly 
speaking, for non half-filled case we never obtain Mott- 
Hubbard insulator in DMFT at all. In Fig. \TE\ we show 
the density of states of Anderson-Hubbard model with 
electron concentration n — 0.8 for different degrees of 
disorder A and U = 6. 01?, representing typical case of 
the doped Mott insulator. Quasiparticle band overlaps 
now with lower Hubbard band and is smeared by dis- 
order, which is just the expected behavior in metallic 
state. Nothing spectacular happens also with conductiv- 
ity, which is shown for the same set of parameters in Fig. 
1161 It shows typical behavior associated with disorder 
induced Anderson MIT. Small signs of transitions to the 
upper Hubbard band can be seen for u> ~ U (see insert 
in Fig. ITS]) . Thus, doped Mott insulator with disorder is 
qualitatively quite similar to the case of disordered cor- 
related metal discussed above. 



V. CONCLUSION 

We used a generalized DMFT+S approach to calcu- 
late basic properties of disordered Hubbard model. The 
main advantage of our method is its ability to provide 
rather simple interpolation scheme between rather well 
understood cases of strongly correlated system (DMFT 
and Mott-Hubbard MIT) and that of strongly disordered 
metal without Hubbard correlations, undergoing An- 
derson MIT. Apparently this interpolation scheme cap- 
tures the main qualitative features of Anderson-Hubbard 
model, such as the general behavior of DOS and dy- 
namic (optical) conductivity. The overall picture of zero- 
temperature phase diagram is also quite reasonable and is 
satisfactory agreement with the results of more elaborate 
numerical work—. Actually, our DMFT+S approach is 
much less time-consuming than more direct numerical 
approaches, such as that of Ref. 1 -, and in fact allows us 
to calculate all basic (measurable) physical characteris- 
tics of Anderson-Hubbard model. 

Main shortcoming of our approach is its neglect of in- 
terference effects of disorder scattering and Hubbard in- 
teraction, which leads to independence of Anderson MIT 
critical disorder A c on interaction U. Importance of in- 
terference effects is known for a long timei^, but its ac- 
count was only partially successful in case of weak cor- 
relations. At the same time, the neglect of these inter- 
ference effects is the major approximation of DMFT+E, 
allowing to derive rather simple and physical interpola- 
tion scheme, allowing the analysis of the limit of strong 
correlations. Attempts to include interference effects in 
our scheme are postponed for future work. 

Another simplification is, of course, our assumption of 
non-magnetic (paramagnetic) ground state of Anderson- 
Hubbard model. The importance of magnetic (spin) ef- 



fects in strongly correlated systems is well known, as well 
as the problem of competition of ground states with dif- 
ferent types of magnetic ordering 8 . Importance of disor- 
der in the studies of interplay of these possible ground 
states is also quite evident. These may also be the sub- 
ject of our future work. 

Despite these shortcomings, our results seem very 
promising, especially concerning the influence of strong 
disorder on Mott-Hubbard MIT and the overall form of 
the phase diagram at zero temperature. The changes 
in phase diagram at finite temperatures will be the sub- 
ject of further studies. Non-trivial predictions of our 
approach, such as the general behavior of dynamic (op- 
tical) conductivity and, especially, the prediction of dis- 
order induced Mott insulator to metal transition can be 
the subject of direct experimental verification. 
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APPENDIX A: EQUATION FOR RELAXATION 
KERNEL 

Let us follow the standard approach of self-consistent 
theory of localizatio n 14 i 15 ' 16 i 17 ' 18 i 19 , taking into ac- 
count the DMFT contributions T, R,A (e) in single-particle 
Green's functions ([30)) and not limiting ourselves to the 
usual limit of small u>. 

Consider Bcthc Salpctcr equation relating the full two- 
particle Green's function 3> p p, (cj, q) to irreducible vertex 
[/pp, (a;, q), accounting only for impurity scattering in 
vertices, but built upon Green's functions given by (|30|) . 
This equation can be written now as a generalized kinetic 
equation of the following for m 14 i 15 i 16 i 17 i 18 i 19 : 

(i}-e(p)-AI^,( W ))$^( W) q) = 

= -ag p +£^(",q)*^(",q)J (Ai) 

where AG P = G R (e+,p+) - G A (e_,p_). The main dif- 
ference with similar equation of Refs i 14 ' 15 i 16 ' 17 i 18 ' 19 is the 
replacement u> — > Co. 
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Let us sum both sides of ()A1 1) and of the same equation 
multiplied by (pq) (where p = and 4 = j^y are appro- 
priate unit vectors) over p and p', taking into account 
an exact Ward identity^: 



^t/^( W) q)AG p - 



(A2) 



and using an approximate representation (cf. Ref. 



14 \ 



(«, q) + 



AG p (pq) 0RA 



EpAG p (pq) 



ri"|2 Is 



l (w,q) (A3) 



where $° flj4 (o;,q) = E PP ' ^ P p^( w J q) is our 1°°P 
I, while <&?f A (^q) = Epp,(M)*^(w,q). Im- 
portant difference from similar representation used in 
Refs^^^JI^^ is that JAU) is not limited to small 



Now (for q — > 0) we obtain the following closed system 
of equations defining both $® RA (ui, q) and $5^(0)^): 

wfc^w, q)- < t> > A (w, q) = - AG p ( A4 ) 



(W + M( W ))$?f A (^, q) - liL^g*°fA(a,, q) = 



where relaxation kernel is given by: 



^ ,pp- > i-i — 1>> P p ^■'-i'-"jp , ( 



V . (pq) AGpEO, q) AG p , (p'q) 



EpAGp 

with average velocity < u > defined as: 



(A5) 



Ep|v P |AGp de(p) 

<u> =i^AGp- ;Vp =V' (A6) 

From (IA4I) we immediately obtain: 



-EpAGp 

u) + iD{uj)q 2 



(A7) 



which for small u> reduces to ([33]) with generalized diffu- 
sion coefficient given by ([38]) . 



Using for irreducible vertex £7p5.^(w,q) an approx- 
imation of "maximally crossed" diagrams and in- 
troducing the standard self-consistency procedure of 
Refs i 14 ^ 15 ' 16 ^ 17 ' 18 ^ 19 (i.e. replacing the Drude diffusion co- 
efficient in Cooperon contribution to irreducible vertex 
by the generalized one defined by ([35])), we obtain from 
(|A5[) our expression (]39p for relaxation kernel. 

Our equation (|4"T]) for the generalized diffusion coeffi- 
cient (which is in general case complex) reduces just to 



the usual transcendent equation, ft was solved by iter- 
ations for every value of Q, taking into account that for 
d = 3 and cutoff given by ([42]) . the sum entering ([4T]) 
reduces to: 



1 



1 



u + iD{uj)q 2 2tt 2 iD{uj)k 2 J y 2 + 



y 2 dy 



« 



2tt 2 iD{u)k 2 



iD{uj)kl 



arctg 



iD(ui)kl 

(A8) 
'iD{uj)kl 



For finite frequencies w we use *° flj4 (q, a)) given by 
(|A7|) , so that an expression ([25]) for dynamic conductivity 
is to be rewritten as: 

Rea(u;) = ^ £ de[f(e.)-f(s + )]x 



AS 



RR, 



The second term here was taken in the "ladder" approx- 
imation: 



E p G«( e+! p + )G^( e _ !P _) 



l_ A 2E p G«( £+ ,p + )G«( £ _ )P _) 

(A10) 

This contribution (non singular at small uj) is irrelevant 
for conductivity at small u — > 0, but leads to finite cor- 
rections with increasing u>. Eq. (|A9[) is our final result, 
which was analyzed numerically in a wide interval of fre- 
quencies (for small uj it reduces to ([37]) ). 



APPENDIX B: "BARE" ELECTRON 
DISPERSION AND VELOCITY 

We consider the "bare" energy band with semi-elliptic 
DOS (|43p. Assuming isotropic electron spectrum e(p) = 
e(|p|) = e(p) and equating the number of states in spher- 
ical layer of momentum space to the number of states 
in an energy interval [e, e + cfe]], we obtain differential 
equation determining the energy dispersion e(j>): 



4irp 2 dp 
(2tt) 3 



N (e)de 



(Bl) 



For quadratic energy dispersion e(p) close to the lower 
band edge we get the initial condition for Eq. (jBip for 
p — > and e — * —D. Then we obtain 



6ir ( 7r — iy5 + — sin(2(p) 



(B2) 



with </? = arccos(4y) and momentum in units of inverse 



lattice spacing. Eq. (]B2p implicitly defines "bare" energy 



dispersion e(p) for electronic part of the spectrum e S 

[-AO]. 

For half-filled band we easily determine the Fermi mo- 
mentum as: 



p F = p(e = 0) = (3tt 2 ) 3 
We also need electron velocity |v p | = 



entering into the expression 
From (IB1I) we obtain: 



9e(p) 



dp 



(B3) 

= ac(£) 
9p 



6|) for average velocity. 



P 



1 



de _ 

rf^ ~~ 2tt 2 jV Q (e) 



(B4) 



where p is given by Eq. (IB2j) . 

To obtain quadratic dispersion for hole part of the 
spectrum (e e [0, D}) close to the upper band edge 
(e — * D) we introduce the hole momentum p = 2pp — p 
and write similarly to (|B 1[) : 



Putting p — > at the upper band edge e — > 0, we get: 



P 



6n (tp — —sin(2ip) 



(B6) 



Then for velocity at the hole part of the spectrum we 
obtain: 



de 



v n = — = — - = 



dp 



de 
dp 



~2 
P 



1 



2TT 2 N (e) 



(B7) 



Eqs. (|B4[) and (|B7|) define energy dependence of |v p |. It 
is easily seen that velocity is even in energy and becomes 
zero at the band edges. These expressions allow us to 
change from momentum summation (e.g. in Eq. (|A6|) ) 
to energy integration. 



4irp 2 dp 
(2tt) 3 



-N (e)de 



(B5) 
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FIG. 3: Bethe-Salpeter equation for polarization loop in 
DMFT+E approach. Circle represents irreducible vertex part 
in particle-hole channel of DMFT, which contains only local 
Hubbard interactions. Unshaded rectangular vertex repre- 
sents corrections from impurity scattering only, implicitly in- 
cluding the case free particle-hole propagation. 



FIG. 1: Typical "skeleton" diagrams for the self-energy in 
the DMFT+S approach. The first two terms are examples of 
DMFT self-energy diagrams; the middle two diagrams show 
contributions due to random impurity scattering represented 
as dashed lines. The last diagram (b) is an example of ne- 
glected diagram leading to interference between the local Hub- 
bard interaction and impurity scattering. 
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FIG. 2: Full polarization loop with vertex part, describing all 
interactions and impurity scatterings in particle-hole chan- 
nel. The loop without vertex corrections is included implic- 
itly. Here p± = p ± § , e± = e ± f . 



FIG. 4: Effective vertex 7 fe (icj,q = 0) used in calculations of 
conductivity. 
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FIG. 7: (Color online) Density of states of half-filled 
Anderson-Hubbard model for different degrees of disorder A 
and U = 4.5D, typical for Mott insulator. 



FIG. 5: Diagrammatic representation of (a>,q) and 




0,0 



FIG. 8: (Color online) Restoration of quasiparticle band by 
disorder in coexistence (hysteresis) region for U = 2.5D, ob- 
tained from Mott insulator with decreasing U. 

FIG. 6: (Color online) Density of states of half-filled 
Anderson-Hubbard model for different degrees of disorder A, 
and U = 2.5D, typical for correlated metal. 




FIG. 9: (Color online) Real part of dynamic conductivity for 
half-filled Anderson-Hubbard model for different degrees of 
disorder A, and U = 2.5D, typical for correlated metal. Lines 
1,2 are for metallic phase, line 3 corresponds to the mobility 
edge (Anderson transition) , lines 4,5 correspond to correlated 
Anderson insulator. 



FIG. 11: (Color online) Real part of dynamic conductivity 
of half-filled Anderson-Hubbard model for different degrees 
of disorder A and U = 1.515, comparison of self-consistent 
theory (full lines) with "ladder" approximation (dotted lines). 
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FIG. 10: (Color online) Real part of dynamic conductivity of 
half-filled Anderson-Hubbard model for different degrees of 
disorder A and U — 4.5D, typical for Mott insulator. Lines 
1,2 correspond to Mott insulator, line 3 corresponds to the 
mobility edge (Anderson transition), lines 4,5 are for corre- 
lated Anderson insulator. Insert - enlarged region of small 
frequencies. 



FIG. 12: (Color online) Real part of dynamic conductivity of 
half-filled Anderson-Hubbard model for different degrees of 
disorder A and U = 2.5D, obtained from Mott insulator with 
decreasing U. 




FIG. 13: (Color online) Zero temperature phase diagram of 
paramagnetic Anderson-Hubbard model. Borders of Mott 
insulator phase (7 c i, C 2(A) are shown as obtained from Eqs. 
(|46p . while different symbols show values calculated either 
from DOS or conductivity behavior. Dotted line defines the 
border of coexistence region obtained with decreasing U from 
Mott insulator phase. Line of Anderson transition is given by 
calculated value of A c = 0.37. 



FIG. 15: (Color online) Density of states of Anderson- 
Hubbard model with electron concentration n — 0.8 for dif- 
ferent degrees of disorder A and U = 6.0D, representing the 
doped Mott insulator. 
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FIG. 14: (Color online) Disorder dependence of static conduc- 
tivity, obtained for several values of U and showing disorder 
induced Mott insulator to metal transition. At the insert we 
show static conductivity dependence on U close to Mott tran- 
sition, including typical hysteresis behavior obtained with U 
decreasing from Mott insulator phase. 



FIG. 16: (Color online) Real part of dynamic conductiv- 
ity of Anderson-Hubbard model with electron concentration 
n = 0.8 for different degrees of disorder A and U = 6.0 D, rep- 
resenting the doped Mott insulator. Insert - high frequency 
behavior with signs of transitions to the upper Hubbard band. 



